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На основе экспериментально определенной зависимости выхода реакции от исходной концентрации активного газа и зависи- 
мости скорости расширения паро-капельного облака предложен подход к анализу кинетики химических реакций при электри- 
ческом взрыве проводника АІ в активных газах. 


Введение 

Электрический взрыв металлической прово- 
лочки (ЭВП) используется для получения как вы- 
сокодисперсных металлов, так и их соединений с 
химически активными газами (0 2 , N 2 , ЫН„ углево- 
дороды и др.) [1-4]. Несмотря на стехиометричес- 
ки достаточные количества этих газов в электровз- 
рывной камере, химические соединения зачастую 
образуются в виде композита химическое соедине- 
ние/металл или даже примесей к металлу. Состав 
композита зависит от свойств реагентов и условий 
эксперимента - введенной энергии, концентрации 


реагентов, давления в системе. Так как состав ком- 
позита важен для формирования свойств материа- 
лов, создаваемых на основе получаемых этим спо- 
собом нанопорошков, то определение причин и ус- 
ловий образования тех или иных составов и чистых 
веществ является актуальной задачей. 

Специфика реакций, происходящих при ЭВП, 
состоит в том, что они начинаются при очень высо- 
кой начальной температуре (~10 4 К) и происходят 
при ее быстром снижении (~Ю 8 К/с), то есть систе- 
ма первоначально находится в состоянии "холод- 
ной" плазмы (ионизация газа ~1 %) с последующей 
эволюцией к обычным условиям. 
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В работах [1,4] были исследованы зависимости 
выхода реакций от концентраций и давлений газов. 
На основании проведенных исследований в рабо- 
тах [4, 5] предложен подход к анализу кинетики и 
механизма реакций, происходящих при ЭВП. В 
предложенном подходе ионизация не учитывается, 
для анализа используется зависимость выхода ре- 
акции от начальной концентрации, скорость рас- 
ширения продуктов взрыва проволочки использу- 
ется для определения температуры реагентов. Со- 
поставление этих данных позволяет найти такие 
параметры системы, как радиус, объем и темпера- 
тура, при которых прекращается реакция. В этих 
работах указано, что более подробное рассмотре- 
ние возможно с использованием двух моделей: 

(а) модель мгновенного перемешивания реагентов, 
при которой перемешивание пара металла с ак- 
тивным газом начинается сразу после начала 
разлета продуктов ЭВП, реакция начинается 
после снижения температуры системы ниже 
верхнего порога реакции и степень превраще- 
ния проволочки определяется объемом газа, 
доступного для перемешивания за время сниже- 
ния температуры до нижнего порога реакции; 

(б) модель расширения паров металла с последую- 
щей реакцией в тонком слое, в которой предпо- 
лагается предварительное концентрирование 
паров металла в цилиндрическом слое [9] с пос- 
ледующим взаимодействием и перемешивани- 
ем, аналогично модели (а). 

Кроме того, возможен третий механизм: 

(в) перемешивание паров металла с активным га- 
зом сначала не происходит, реакция начинается 
на более поздних стадиях расширения продук- 
тов ЭВП, когда весь металл сконденсировался в 
капле. Реакция идет первоначально в монослое 
на поверхности капель металла, затем, вслед- 
ствие саморазогрева, более мелкие капли сгора- 
ют, а крупные - лишь покрываются тонким сло- 
ем продукта, так как их саморазогрев меньше; в 
этом случае степень превращения металла оп- 
ределяется распределением по размерам пер- 
вичных металлических частиц. 

В настоящей работе проведен анализ экспери- 
ментальных данных по зависимости выхода реак- 
ций алюминия с азотом, аммиаком и кислородом с 
использованием первой модели. Анализ моделей 
(б) и (в) будет проведен в последующих работах. 

1. Экспериментальные данные 

Экспериментальные зависимости выходов хи- 
мических реакций (77) от начальных концентраций 

активных газов (с 0 ) были получены на установке, 
описанной ранее [1]. Электрический взрыв алюми- 
ниевой проволоки диаметром с/=0,38 мм и длиной 

/=100 мм, массой 0,03051 г (1,13 ммоль) производи- 
ли в атмосфере, содержащей активные газы (0 2 , 

ЫН 3 , N 2 ) при давлениях 1...4-10 1 * * * 5 Па. При этом ввод 
энергии І?=1,8...2,0Д. (Е с - энергия сублимации) 


производили за время порядка 2 мкс. Методики 
анализа продуктов описаны в [5]. Результаты экс- 
периментов приведены на рис. 1. 

Концентрация N2, % 



Рис. 1. Зависимость выходов реакций АІ с Оі+Аг (кривая 1), 
N Н з +Аг ( кривая 2) и Ыі+ Ат ( кривая 3) от р активных газов 

Более подробно была исследована система 
А1+1Ч,. Для нее, в частности, были проведены иссле- 
дования продукта реакции А1+]Ч 2 методом элект- 
ронной микроскопии высокого разрешения 
(НКЕМ). При применении дифракции электронов 
и метода светлого/темного поля было показано, что 
в частично азотированных образцах (®(А1М)«70 %), 
металлические ядра отсутствуют, а на поверхности 
частиц АШ имеются блоки типа дендритов той же 
структуры, что и основная частица. Высокое разре- 
шение также не позволило обнаружить иную решет- 
ку, то есть неокисленный металл, видимо, находит- 
ся в аморфном состоянии или в виде отдельных 
крупных частиц, не попавших в поле микроскопа. 
(В отличие от этих данных, для других систем часто 
наблюдаются продукты в виде структур ядро-обо- 
лочка [5]). 

2. Модель реакции 

С использованием наших и литературных дан- 
ных [6-10] для систем, в которых реакция идет в га- 
зовой фазе, может быть предположена следующая 
картина. 

- После нарушения сплошности проводника 
происходит расширение, диспергирование и 
испарение металла с образованием паро-ка- 
пельного облака [6] квазицилиндрической фор- 
мы [9, 10]. Будем далее рассматривать вариант 
процесса, в котором существует квазиравнове- 
сие между паром и каплями металла и расход 
пара на реакцию приводит к немедленному ис- 
парению металла из капель, так что лимитирует 
скорость процесса химическая стадия, а стадия 
испарения-конденсации является быстрой. 

- Предположим далее, что после снижения тем- 
пературы системы до температуры, близкой к 
Гран, (при К ѵш =\ (химической реакции (но нес- 
колько выше ее), происходит образование про- 
дукта и, вследствие этого, быстрое перемешива- 
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ние пара и активного газа, так что объемную на- 
чальную концентрацию активного газа перед 
началом реакции можно принять равной его ис- 
ходной концентрации (с 0 ). 

- Будем рассматривать только газо-фазную тер- 
мически активируемую реакцию между метал- 
лом и газом, так как ее скорость на много поряд- 
ков больше, чем реакция на границе конденси- 
рованных фаз вследствие ее торможения диф- 
фузией реагента через слой продукта; это поло- 
жение соответствует экспериментальным дан- 
ным для системы А1+М 2 , указывающим на отсу- 
тствие частиц типа ядро-оболочка, которые об- 
разовались бы в случае гетерогенной реакции. 

- Предполагаем, что реакция прекращается в ре- 
зультате понижения температуры до некоторого 
порогового значения температуры Т„, которая 
ниже 7^ авн ; этой температуре соответствуют т„, г„ 
и Ѵ„ - время, радиус и объем цилиндра, при ко- 
торых прекращается реакция. Для гетерогенной 
реакции такие параметры тоже могут быть най- 
дены, однако они будут иметь смысл эффектив- 
ных параметров. 

3. Расчет распределения температур 
в зоне химической реакции 

Прямое измерение температуры в процессе 
ЭВП представляет значительную трудность. В ос- 
новном при измерениях высоких температур ис- 
пользуется метод пирометрии. Результаты измере- 
ний этим методом показывают, что при ЭВП раз- 
виваются температуры свыше ІО 4 К [6]. Однако 
пирометрия, видимо, дает завышенные значения 
температуры для атомов и молекул вследствие тер- 
мической неравновесности процесса на первой его 
стадии [11]. В этом случае понятие температуры к 
системе в целом, как известно, неприменимо; и 
этим параметром можно пользоваться лишь с оп- 
ределенными оговорками - обычно говорят об эф- 
фективной температуре. 

Если все же предположить, что вводимая энер- 
гия (2Е&605 кДж/моль, с = 21,39 Дж/(моль-К) -для 
А1 [12]) целиком расходуется на нагревание и испа- 
рение всего металла, то температура системы в слу- 
чае А1 будет порядка 14000 К. В действительности 
полного испарения металла не происходит, а обра- 
зуется паро-капельная смесь. В результате этого 
часть энергии, направленная на испарение метал- 
ла, расходуется на нагрев пара, и реальная темпера- 
тура системы может быть выше. С другой стороны, 
при ЭВП имеет место генерация ударных волн, на 
которую может уходить до 50 % введенной в про- 
водник энергии [13], кроме того, часть энергии те- 
ряется на излучение. Таким образом, оценка тем- 
пературы в зоне ЭВП представляет сложную задачу. 
Необходимо также учитывать, что происходящие 
при ЭВП химические реакции выделяют или пог- 
лощают энергию, что также изменяет температуру. 

Наиболее простым способом эффективную 
температуру в зоне ЭВП можно рассчитать из ско- 


рости движения контактной границы паро-капель- 
ное облако - окружающий газ, в предположении, 
что передаваемая энергия первоначально расходу- 
ется на разрыв химических связей и одностороннее 
движение атомов. При этом можно воспользовать- 
ся законом Максвелла [14]: 


Т. = 


к пгіѴ 
8Я 


( 1 ) 


где Т э - эффективная температура одномерного 
движения молекул, т - молекулярная масса пара, 
]Ѵ - скорость движения контактной границы, 
К - газовая постоянная. 


Изменения средних значений молекулярных 
масс и квадратов скоростей молекул в формуле (1) 
для разных систем при постоянстве введенной энер- 
гии компенсируют друг друга. Таким образом, Т э за- 
висит только от введённой энергии и будет одинако- 
ва для всех систем. Температура системы после пере- 
мешивания и установления равновесия между раз- 
личными степенями свободы будет равна Т,„ делен- 
ной на общее число степеней свободы и расчетная Т 
системы будет зависеть от числа степеней свободы. 

Для систем 4А1+30 2 и 2А1+1Ч 2 суммарные числа 
степеней свободы составили 4,98 и 4,26 соответ- 
ственно. В случае аммиака суммарное число степе- 
ней свободы было рассчитано для системы 
А1+1/21Ч 2 +3/2Н 2 (т.к. при температуре выше 
1000... 1200 К весь аммиак находится в диссоцииро- 
ванном состоянии) и составило 5,18. 


Второй метод основывается на определении из- 
менения температуры газа в результате адиабати- 
ческого сжатия. Адиабатическое сжатие газа, про- 
исходящее при прохождении ударной волны, при- 
водит к увеличению температуры газа, которое 
можно оценить по уравнению [15]: 


Т 


27 ( 7 - 1 ) 

(7+1) 2 


м- 


( 2 ) 


где Т 0 - начальная температура газа (300 К), 
М- число Маха ( Ѵ/Ѵ щш ), у=с р /с ѵ - показатель адиа- 
баты. Известно [9], что при ЭВП сначала образует- 
ся и проходит ударная волна, а затем идет расшире- 
ние паро-капельной смеси, которое также вызыва- 
ет адиабатическое сжатие. Температура может быть 
рассчитана как по отношению к газу в ударной вол- 
не, которая опережает фронт паро-капельного об- 
лака, так и для самого этого фронта - в зависимос- 
ти от используемых данных для скорости. В данной 
работе расчеты проведены для фронта расширяю- 
щегося облака атомов А1 и для смесей его с газами. 

Достаточно подробные данные по скорости 
разлета продуктов ЭВП в воздушной атмосфере 
приведены в [1]. Зависимости радиуса разлета для 
алюминия в воздухе и аргоне приведены в [16]. Из 
работы [16] видно, что скорость разлета малочув- 
ствительна к природе газа при одинаковых Е/Е с . 

Для расчета зависимости температуры от време- 
ни разлета использована экспериментальная кри- 
вая скорости разлета продуктов при ЭВП [1], рис. 2, 
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кривая 1. По этой кривой рассчитали зависимость 
Т 3 от времени (г) по формуле (1), из неё получили 
кривую 2 (А1) для паров А1, а затем рассчитали зави- 
симости Т от т для других систем с учетом среднего 
числа степеней свободы - на рис. 2 , кривые 2 (О,), 
2 (N 2 ) и 2 (МН 3 ), соответственно. Кривая 3 на этом же 
рисунке - результат расчета по ур. (2) для паров А1. 

Из приведенных температурных зависимостей 
видно, что температурная кривая 3, полученная по 
ур. (2), лежит близко к кривой 2(А1) для паров А1, 
найденной по ур. ( 1 ). Для дальнейших расчетов ис- 
пользованы кривые, полученные по ур. ( 1 ), т.к. в 
этом случае учтено различие состава систем. 

4. Оценка параметров реакционного цилиндра 
по значениям концентраций активных газов, 
необходимых для полного превращения металла 

В соответствии с принятой моделью реакция 
прекращается в результате снижения температуры 
в некотором объеме вокруг диспергирующегося 
проводника. Зная температурный профиль Т(т), 
выход реакции ( 77 ), количество вступившего в ре- 
акцию газа (ѵ), а также длину и массу проводника, 
пользуясь данными рис. 1 и 2 и формулами: 

ѵ=Р Ѵп/КТ, Ѵ=я-г}1, (3) 

можно вычислить параметры условного реакцион- 
ного цилиндра, такие как объем Ѵ„, время т„, темпе- 
ратуру Т„ и радиус прекращения реакции г„. Рассчи- 
танное из длины и диаметра проводника для наших 
условий количество А1 составляет 1,13 ммоль. В таб- 
лице приведены стехиометрические количества 
вступающих в реакцию газов (ѵ 0 для одного взры- 
ва). В строках 1-3 приведены результаты расчетов с 
использованием уравнения ( 1 ) размера цилиндра, 
температуры и времени прекращения реакции в 
предположении, что в реакционном цилиндре вы- 
горает весь газ. Давление, при котором происходит 
100 %-ное превращение металла в этих газах, опре- 
делено из рис. 1 и по данным [ 1 ]. 



1 2 3 4 5 6 х > мкс 


Рис. 2. Зависимость скорости движения контактной границы 
паро-капельное облако — окружающий газ и темпе- 
ратуры от времени: 1) скорость движения контактной 
границы [1]; 2) зависимости температуры от времени 
для систем, расчет по ур. (1): А1+0 2 (2(О г )), А!+М г 
(2( N 2 )), АІ+МН з (2(ЫН 3 )), АІ (пар) (2(АІ)), 3) зависи- 
мость температуры от времени, расчет по ур. (2); 
4) радиус (К) контактной границы паро-капельное 
облако — окружающий газ [1] 


Из этих данных видно, что радиус и температу- 
ра прекращения реакции, рассчитанные для азота, 
физически не корректны: реакция по этому расче- 
ту прекращается уже на участке разгона паро-ка- 
пельной смеси, в области роста температуры. Сама 
Т„ выше температуры равновесия диссоциации АПЧ 
(3795 К [5]). Это объясняется грубостью принятой 
оценки, которая не учитывает кинетики превращения. 

Таблица. Значения параметров для различных систем, оп- 
ределенные из экспериментальных зависимостей 
выхода от концентрации активного газа (при 
г]=Ю0 % и р=((с 0 ), Рио — парциальное давление 
активного газа, необходимое для полного прев- 
ращения металла) 


Система и 
условия 
расчета 

К' 

ммоль 

І'ох.мл 

К, мл 

С Ю" Па 

т п , МКС 

/-„•ІО 2 , м 

Т п , К 

АІ+0 2 +Аг; 
/7=100 % 

0,85 

95,2 

95,2 

0,2 

>6,2 

1,72 

зоо 

АІ+Щ+Аг; 

/7=100% 

1,13 

12,7 

12,7 

2 

3,6 

0,62 

2040 

АІ+і\і 2 
/7=100 % 

0,567 

0,846 

0,846 

15 

2,1 

0,15 

4870 

АІ+і\і 2 +Аг; 
ц =Ас) 

0,567 

ѵ ш =Ас) 

5,60 

15 

2,84 

0,404 

3710 


Однако приведенные значения г„ и Т„, рассчитан- 
ные из экспериментальных данных по полному выходу 
химических соединений, в сравнительном плане впол- 
не соответствуют обычным представлениям о хими- 
ческой активности газов по отношению к алюминию. 
Эта активность уменьшается в ряду О,— >МН 3 ->М 2 по 
совокупности причин, связанных с термодинамикой 
реакций и строением молекул реагента: О, - биради- 
кал, из-за наличия неспаренных электронов он имеет 
повышенную реакционную способность и наиболь- 
шее сродство к алюминию; N 2 - очень прочная моле- 
кула с тройной связью, известная своей кинетической 
инертностью. Молекула 1ЧН 3 по совокупности этих 
признаков должна занимать промежуточное положе- 
ние. Это совпадение рядов - простое отражение того, 
что выход реакции возрастает с увеличением ее ско- 
рости в принятой геометрии реакционной зоны. 

5. Оценка параметров реакционного цилиндра 
с использованием кинетических уравнений 
газофазных реакций 

Поскольку размер реакционного цилиндра в 
принятой модели определяется температурой 
прекращения реакции и реакция происходит в од- 
ном и том же температурном интервале независи- 
мо от концентрации газа, то скорость газовой реак- 
ции можно записать в обычном виде: 

-^ = кс", (4) 

а т 

где с - концентрация газа (концентрация паров АІ 
исключена, так как она не зависит от степени прев- 
ращения вследствие равновесности стадии испаре- 
ния-конденсации); к(Т) - зависящая от температу- 
ры константа скорости реакции. 
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Технические науки 


Выход реакции равен: 


'7 = 



(5) 


где ѵ 0 и ѵ - количество газа, стехиометрически не- 
обходимое и оставшееся после реакции. Т.к. ѵ=с V, 
то можно выразить через концентрацию, объем га- 
за, необходимые для полного превращения про- 
водника (Ѵ сш ), и Ѵ„. При этом величина Ѵ„ постоян- 
на для каждой химической системы при неизмен- 
ных условиях взрыва. Рассмотрим два случая: 


при п = 2 : 


п = - 


ѵ„с: 


г 


Ѵ 1 + С о^У 


( 8 ) 


Зависимости экспериментальных значений г](с 0 ) 
для различных порядков реакции А1 с N 3 в линеари- 
зующих координатах представлены на рис. 3, а-в. 

Видно, что наилучшее согласие эксперимен- 
тальных и теоретических зависимостей имеет мес- 
то при / 7 = 0 , 5 . 


1. "Медленное" охлаждение - Ѵ„> Ѵ т „ тогда 

с 'V с 

ц = 1 где с = 9 Р 

Ѵ “ {\+С'(п-\)Р\-' 


- для шП и с=с 0 ехр(-Т) - для п= 1 ; так что 


'7 = 



- для / 7 = 0,5 и 

<7 = 1 -—с 0 е Р 

«о 

- для п= 1 . 

Т„ 

Здесь Р = | к(т)сІт - не зависящая от концентра- 

О 

ции функция, характеризующая реакционную спо- 
собность системы, где к(т) - изменяющаяся во вре- 
мени (с изменением температуры) константа ско- 
рости реакции. 

Экспериментальные данные (рис. 1) показыва- 
ют, что этот случай не реализуется, теоретическая 
зависимость ( с 0 ) противоположна эксперименталь- 
ной. 


2. "Быстрое" охлаждение - Ѵ„< Ѵ сш , тогда 
V. 


Ч = — (с 0 -с), и 
ѵ„ 


при и = 0,5: 


при п = 1 : 


77 = 


Ѵ*Г< РУ, 

и. 4и„ 


V г 

ц =— (1-е-О; 


( 6 ) 


( 7 ) 


Графический анализ дал следующие значения: 
/=0,308; К„=5,6 мл. Рассчитанные также по ур. (3) 
и рис. 2 величины г„, Т„ и т„ приведены в таблице 
(строка 4). В отличие от данных, приведенных в 
строках 1-3, этот расчет дал физически корректные 
данные: реакция прекращается на участке пониже- 
ния температуры. 

При этом порядок реакции и=0,5 соответствует 
предполагаемому механизму: 

N^2 К,К=к 1 /к;=[ЩУ[Щ-, (А) 

N + АіДаііМ, (Б) 

где к ь к[, к 2 - константы скоростей соответствую- 
щих прямых и обратной реакций, К - константа 
равновесия диссоциации азота. 

В соответствии с этим механизмом стадия (А) 
обратима, а (Б) - лимитирующая, причем концент- 
рация [А1] не зависит от реакции, поэтому скорость 
пропорциональна корню квадратному из концент- 
рации N 3 : 

^ АІЫ] = к 2 ■ [А1] ■ [N1 * к 2 ■ [А1] К 0,5 ■ [Ы ,] °’ 5 . 
сіі 

Другой важный параметр - температура прек- 
ращения реакции Т„ в этом случае меньше темпера- 
туры равновесия реакции (Б). Это делает вероят- 
ным газофазный механизм, но не исключает пол- 
ностью другие, т.к. приведённые оценки темпера- 
турных профилей нельзя считать точными, не ясен 
также механизм перемешивания газов за времена 
порядка 3 мкс и на расстояниях около 4 мм. Эти 
вопросы нельзя решить на основании представлен- 
ных данных и расчетов. 

Выводы 

Таким образом, если химическая реакция проте- 
кает в газовой фазе, то предложенная методика поз- 
воляет на основании экспериментальных данных о 
выходах химических реакций и скоростях расшире- 
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Рис. 3. Проверка соответствия экспериментальных данных для системы АІ+И 2 кинетическими ур. ( 6—8) 
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ния паро-капельного облака оценить такие пара- 
метры, как температуру, радиус, объем и время 
прекращения реакции, а также функцию реакцион- 
ной способности Р(к,т) и порядки реакций. Эти па- 
раметры, в свою очередь, могут быть использованы 
для обсуждения соответствия предложенному меха- 
низму с точки зрения реалистичности (корректнос- 
ти) их самих, или следующих из них выводов. Выше 
уже отмечена качественная корректность всех най- 
денных параметров, кроме величины Г, которая не 
обсуждалась. Эта величина имеет смысл интеграль- 
ной константы скорости реакции. Она включает в 
себя независящие от температуры и времени вели- 
чины энергии активации и предэкспоненциального 
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